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We extend an approximation earlier developed by us for the single-impurity Anderson model to a 
full-size impurity solver for models of interacting electrons with multiple orbitals. The approximation 
is based on parquet equations simplified by separating small and large energy fluctuations justified 
in the critical region of a pole in the two-particle vertex. We show that an Z-orbital model with 
most general interaction is described within this approximation by 4l 2 x 4l 2 matrices and is Fermi 
liquid in the metallic phase. We explicitly calculate properties of a paramagnetic solution of a two- 
orbital Hubbard model with a Hund exchange and orbital splitting within the dynamical mean-field 
approximation. We trace the genesis of a metal-insulator transition induced by a crystal field and 
vanishing of the Kondo quasiparticle peak in strongly correlated orbitally asymmetric systems. 
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I. INTRODUCTION 

Understanding the effects of electron correlations is 
important not only for comprehending the behavior of 
a wide class of materials. It also allows one to get to 
fundamental principles of quantum theory. Quantum dy- 
namical effects caused by electron correlations are still in 
the center of interest of condensed-matter physicists. Al- 
though a lot of aspects of electron correlations have been 
successfully explained, there are phenomena the picture 
of which is incomplete or lacks substantial pieces. The 
main obstacle in reaching a complete knowledge of the ef- 
fects of electron correlations is the difficulty of a reliable 
description of dynamical many-body phenomena. 

The major effort in grasping strong electron corre- 
lations is focused on single-impurity Anderson model 
(SIAM). 1 This model was conceived to provide a sim- 
ple picture of formation of local magnetic moments. It 
appeared rather soon that formation of local moments is 
connected with the Kondo effect^ Since then SIAM has 
been attracting standing attention of condensed-matter 
theorists and experimentalists. The first wave of interest 
in SIAM was crowned by finding an exact ground state 
of the model with a flat band of conduction electrons for 
all interaction strengths^. The Kondo effect in SIAM was 
then explained and quantitatively described. The exact 
expression for the Kondo asymptotics in SIAM has be- 
come a hallmark and reference for dynamical effects of 
strong electron correlations. 

Impurity models suit for a description of quantum 
dots and nano-particles^^ but are unable to describe 
bulk, macroscopic effects of electron correlations. It is, 
however, generally believed that a local electron interac- 
tion in the tight-binding scheme offers enough space for 
the description of most of the dynamical effects of elec- 
tron correlations. Success of the dynamical mean-field 
theory (DMFT) of electron correlations has proved this. 6 
Dynamical mean-field approach, that is local dynamical 
approximations, renewed interest in SIAM, since the ba- 
sic ingredient of DMFT is an impurity solver. The latter 



essentially is an impurity model with a self-consistent 
condition on the local one-electron propagators. The ex- 
act thermodynamic solution of SIAM, however, cannot 
be used, since it does not produce the necessary dynami- 
cal renormalization of the energy spectrum. A new wave 
of attempts to find approximate dynamical solutions of 
SIAM in the strong-coupling regime arose. 

Presently, the most comprehensive quantitative ap- 
proaches producing dynamical properties of the impu- 
rity models are quantum Monte Carlo (QMC) 7 - with its 
continuous-time extension;^ exact diagonalization (ED) 
with a discretized bathji - and the numerical renormal- 
ization group (NRG) . 11 ! 12 The first two approaches aim 
at a precise thermodynamics of the impurity problem, 
while the last one best reproduces low-temperature and 
low-energy scales of impurity models. The first two 
approaches work naturally in the Matsubara formalism 
and need a numerical analytic continuation to determine 
spectral properties, while NRG can be straightforwardly 
extended to dynamical and spectral quantities ^ 3 ' 14 The 
numerical renormalization group represents now-a-days 
the most accurate quantitative low-temperature and low- 
energy dynamical solution of impurity problems. 15 

Apart from pure numerical methods a series of ana- 
lytically controlled approximations have been proposed. 
An approximation based on low-frequency renormaliza- 
tions of Fermi-liquid parameters is able to reproduce 
the Kondo scale in SIAM and fits for heavy Fermion 
systems J£ It is not, however, suitable for the description 
of critical instabilities of the Fermi liquid state. Strong- 
coupling expansions based on the infinite-interaction 
mode l 17 ' 18 reproduce as well the Kondo scale but fail to 
reproduce the Fermi liquid regime in the weak-coupling 
limit. 

Standard truncated weak-coupling expansions fail to 
reproduce correctly the strong-coupling regime. An 
improvement can be reached if the expansion is ap- 
plied on vertex functions. Expansions for two-particle 
vertices are mostly based on multiple scatterings of 
electrons on electrons (T-matrix, TMA) or electrons 
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on holes (random phase, RPA). The non-renormalized 
multiple electron-hole scatterings lead to an unphys- 
ical pole that must be removed by a self-consistent 
treatment. Standardly introduced one-electron self- 
consistency (fluctuation-exchange approximation) re- 
moves the pole but fails to reproduce the Kondo scale* - 9 
A static one-electron self-consistency of the Hartree type 
was introduced in an effective approximate scheme, called 
local-moment approach^ The expansion for the polar- 
ization operator (RPA vertex) uses one-electron propaga- 
tors with a broken spin symmetry in this approximation. 
The spin symmetry of the equilibrium thermodynamic 
state is then recovered via a construction of the self- 
energy calculated from the polarization operator. In this 
way the RPA pole is not crossed and the strong-coupling 
limit keeps the polarization operator in the critical region 
of the RPA pole. The spectral function displays the cor- 
rect Kondo asymptotics of the quasiparticle peak. This 
theory, however, contains an artificial symmetry break- 
ing and an auxiliary fitting parameter, the length of the 
local magnetic moment, in the strong-coupling regime. 

An alternative way of introducing a self-consistency at 
the two-particle level was proposed by us i 21 ' 22 We sug- 
gested to use parquet equations constructed from mul- 
tiple electron-electron and electron-hole scatterings and 
use the criticality of the RPA pole for a simplification 
leading to a manageable theory. This approximation al- 
lows us to control analytically the RPA pole and the gen- 
esis of the Kondo asymptotics. The derived Kondo be- 
havior reproduces universal features of the exact Bethe- 
ansatz solution. 

Real materials are more complex and cannot be de- 
scribed by elementary impurity models. One has to ex- 
tend the methods to multi-orbital Hamiltonians. Not all 
the methods used as impurity solvers for single-orbital 
models work also in the multi-orbital case. The exact so- 
lution can be extended to SU(n) models only in the limit 
U — > ooi 2 ^ There are natural extensions of QMC 2 ^— , 
ED,2L28 and NRG^o methods to multiple-orbitals. All 
the numerically exact approaches are, however, computa- 
tionally rather expensive for realistic atomic structures. 
There is hence still a need for simpler, computationally 
less expensive semi-analytic approaches that could qual- 
itatively correctly describe the transition from the weak 
to the strong coupling regimes in multi-orbital impurity 
and lattice models. 

Some of the existing analytically controlled ap- 
proaches to strong electron correlations have recently 
been extended to multi-orbital models, mostly SU(n) 
Anderson models. Various spectral and thermody- 
namic quantities have recently been calculated within 
the renormalized Fermi-liquid expansion, 31 fluctuation 
exchange^ or local-moment approach* 33 ' 34 The aim of 
this paper is to present a generalization of a two-particle 
approach based on simplified parquet equations to multi- 
orbital models and create a full-scale impurity solver for 
realistic models of correlated electrons. Unlike the most 
of the other multi-orbital approaches, we formulate the 



method so that it can directly be applied not only to 
impurity models but also to lattice models within the 
dynamical mean-field theory. The approximation is for- 
mulated in real frequencies and hence the spectral prop- 
erties are directly available. The low and high frequency 
behavior of the one-electron spectral function, its three- 
peak structure, is qualitatively well reproduced and the 
method allows one to control analytically the critical be- 
havior of Bcthe-Salpeter equations for the two-particle 
vertex. 

The paper is organized as follows. We introduce 
a general multi-orbital perturbation expansion for two- 
particle vertices in Sec. [II] In Sec. IIIII we derive simpli- 
fied parquet equations with electron-hole and electron- 
electron multiple scatterings and show how the self- 
energy is calculated from the full two-particle vertex. We 
present explicit equations to be solved in a two-orbital 
case in Sec. IIVI The derived approximation is used to 
obtain a numerical solution of a two-orbital model with 
Hund's coupling and orbital splitting in Sec. [V] We also 
compare our results with those from other approaches 
there. Reliability and domains of applicability of the de- 
rived approximation are discussed in the final Section lVTl 

II. PERTURBATION EXPANSION FOR 
MULTI-ORBITAL MODELS OF CORRELATED 
ELECTRONS 

We start with a model Hamiltonian with multiple or- 
bitals, where the band structure is resolved and kinetic 
energy is diagonal in orbital indices. The unperturbed 
Hamiltonian then reads 

#o=££e t (k)4(k)c tCT (k) (1) 

k I, a 

where k, a, and i are momentum, spin and orbital in- 
dex of valence electrons. The interaction multi-orbital 
Hamiltonian can generally be represented in the direct 
space as 

^ = lE E ^ 1 2 (xi-x 2 ) C t CTi (x 1 )c.L 2 (x 2 ) 

XI, x 2 ! 
<7l - (T2 1<V,K 

x c KCT2 (x 2 )c A(Ti (xi) (2) 

where all possible inter-orbital transitions are included. 
Orbital-momentum conservation poses a restriction on 
the selection of orbital indices of the interaction matrix 
being l + v = n + A. We use a tight-binding approxi- 
mation taking the coordinate x from the lattice space. 
We will assume the interaction to be local in the lattice 
space, that is, Hubbard-like. A graphical representation 
of the interaction matrix is plotted in Fig. [T] where we 
used latin indices for all the dynamical degrees of free- 
dom characterizing the electron states. They are gen- 
erally four-momenta or Matsubara frequencies in local 
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FIG. 1. General electron-electron interaction in the notation 
we use in this paper. Index a stands for spin, while i,j,k,l 
comprise all other dynamical degrees of freedom character- 
izing the electron state. Conservation law in the dynamical 
variables i + j — k + I holds. 
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FIG. 2. Connecting interactions via electron-electron propa- 
gation. 
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approximations (dynamical mean-field theory) together 
with orbital indices. We use this abbreviated notation 
to simplify summations over intermediate states in the 
multi-orbital perturbation expansion. 

It is necessary to classify all possible interconnections 
of interaction vertices via fermion lines for building up 
the complete multi-orbital perturbation expansion. It is 
not an easy task and we will sum only selected classes of 
two-particle diagrams. We use the graphical representa- 
tion of the interaction vertex from Fig. Q] and choose the 
upper part of it so that spin (not necessarily orbital in- 
dex) is conserved along the horizontal fermion lines. We 
choose electron propagation from left to right and hole 
propagation from right to left. We have three possibilities 
how to connect two interaction vertices by two fermionic 
lines. They correspond to three two-particle scattering 
channels^ 

The first possibility to connect interaction vertices 
is via propagation of two electrons, that is, by parallel 
fermionic lines, cf. Fig. [5J Mathematically this process 
can be represented via a matrix multiplication 

(^G^J^^f&C^^ (3) 

i' 

where the latin indices stand for all dynamical variables 
of the electron (momentum, frequency and orbital index). 
They obey a conservation law j' = k + I — %' . 

Second connection of interaction vertices is via prop- 
agation of an electron-hole pair, that is, by antiparallel 
fermionic lines. Here we must distinguish two options: 
electron-hole scatterings and polarization bubbles. The 
former process is plotted in Fig. [3] and is represented via 
another matrix multiplication 

{UGG.VtfX = £ U^QG^G^U^ (4) 

i' 

where again the orbital-monentum conservation law k' — 
i' + k-l holds. 

A polarization bubble, shielding the interaction, is 



FIG. 3. Connecting interactions via electron-hole propaga- 
tion. 

represented via a matrix multiplication 

a' j' 

(5) 

and is plotted in Fig. |4] Notice that spin in this last 
process is not conserved, but orbital-momentum conser- 
vation restricts the intermediate indices of the dynamical 
variables to k' = j' + i — I. 

We use these elementary scattering processes to con- 
struct the full two-particle vertex V from a two-particle 
perturbation expansion. Once we have this vertex we 
construct the one-electron self-energy as the fundamental 
quantity from which we derive all thermodynamic prop- 
erties. The self-energy is calculated from the two-particle 
vertex via the Schwinger-Dyson equation, diagrammati- 
cally represented in Fig. [5j Mathematically it reads 

j 3 
a' j',k' 

Again the conservation law for intermediate variables 
I' = i + k' — j' holds. The one-electron propagators 
in the perturbation expansion can either be chosen as 
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FIG. 4. Connecting interactions via electron-hole triplet po- 
larization bubbles (vertical channel). 
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FIG. 5. Schwinger-Dyson equation connecting the self-energy 
E with vertex V. The first and the second term on the right- 
hand side of this equation are the Hartree and the Fock con- 
tribution. 



fully dynamically renormalized as suggested by Baym 
and Kadanoff, 3 - or only statically renormalized via the 
Hartree approximation securing that the particle densi- 
ties do not change during the summation of diagrams. 
We discussed earlier that the latter option is more ap- 
propriate in the strong-coupling regime, since we can 
better control the two-particle singularity in the electron- 
hole correlation function and consequently both low and 
high energy scales are qualitatively better reproduced 
than in approximate schemes with the fully renormalized 



propagators. 



III. SIMPLIFIED PARQUET EQUATIONS: 
MATRIX FORMULATION 



The generic function for our perturbation expansion 
is the two-particle vertex T. Its approximate expression 
together with the Schwinger-Dyson equation ([6]) define 
our approximate scheme. We hence start with simple 
approximations on vertex T. We resort our analysis only 
to local approximations, that is, we stay within the dy- 
namical mean-field theory with only local fluctuations 
and Matsubara frequencies, orbital index and spin as dy- 
namical variables. We will treat these degrees of freedom 
separately. 

The simplest approximations for vertex T, beyond 
the bare interaction, are multiple scatterings of the same 
form as defined by elementary processes from the preced- 
ing section. Sums of all multiple scatterings of the same 
type are solutions to Bethe-Salpeter equations. They can 
be represented with the appropriate matrix multiplica- 
tions. For electron-hole and electron-electron multiple 
scatterings we obtain 

$O m ) = U- Ux(iv m ) • ${iv m ) , (7) 
§{iv m ) = U- Lhjj{iv m ) * $(iv m ) , (8) 

respectively. The kernels of these equations are com- 
posed from the bare interaction and either the electron- 
hole bubble or the two-electron propagator. The kernel 
in the electron-hole channel explicitly reads 



with the electron-hole bubble 

Xkcts {Wm) = 4 ^2 G L<Jl {iv m + iu n )G K<72 (mj„) . (10) 

71 

Analogously we construct the kernel of the Bethe- 
Salpeter equation (jHJ with multiple electron-electron 
scatterings. The two-electron propagator then is 

= ^,^2G LrTl {iv' m - iuj n )G vr72 (iuj n ) . (11) 

n 

We denote fermionic and bosonic Matsubara frequencies 
Lj n = (2n + l)nT and v m = 2rmrT, respectively (Boltz- 
mann constant ks — 1). Notice that the transfer fre- 
quency v m , conserved in electron- hole scatterings, is dif- 
ferent from that conserved in the electron-electron scat- 
tering channel, v' rn . If the incoming and outgoing fre- 
quencies of one electron are w n and uj n i , respectively, we 
have v' m = v m + oj n - uv . 

We do not take into account the third channel with 
triplet polarization bubbles. There are two reasons why 
we can afford to do this. First, the effects of the scatter- 
ings in the vertical channel are similar to those from the 
electron-hole channel, at least qualitatively. We hence do 
not lose track of any relevant effect in the strong-coupling 
regime. Second, taking into account only two scattering 
channels simplifies significantly approximate schemes for 
the two-particle vertex.— 

From the solutions of the Bethe-Salpeter equations 
for $ and ^ we obtain the full two-particle vertex 

T(iv m , iv m + ioj n - iw n >) 

= $(ii/ m ) + *(ii/ m + iuj n - iuv) - U . (12) 

Such a vertex defines the self-energy of the so-called 
FLEX approximation if the one-electron propagators in 
the Bethe-Salpeter equations (J7) and © are fully renor- 
malized, that is, they contain the self-energy determined 
by the approximate vertex V from Eq. (TT21 . Such an 
approximation is known to be unreliable in the strong- 
coupling regime^ 

Intermediate electron coupling is marked by proxim- 
ity of a pole in the electron- hole vertex If the transfer 
frequency v = then we can define a small dimensionless 
scale 
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approaching zero in the random-phase approximation 
with T = a solution of Eq. 0, when U -t U c . We 

used Sp A to denote the set of eigenvalues of matrix A. 
This pole is unphysical and indicates that the perturba- 
tion expansion for U > U c must be renormalized. The 
critical interaction U c from RPA defines an upper limit 
on the weak-coupling regime where perturbation theory 
with selected classes of generic diagrams can safely be ap- 
plied. We speak about strong coupling if U > U c , where 
only self-consistent approximations may be reliable. 
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One, and most often used option how to suppress 
the unphysical pole, or better to shift it to infinite in- 
teraction strength is to employ one-particle dynamical 
self-consistency. We use the self-energy calculated from 
the Schwinger-Dyson equation ([6]) with vertex T from 
Eq. (|12[) to renormalize the one-electron propagators in 
Bethe-Salpeter equations and (jHJ. This one-particle 
self-consistency succeeds in shifting the unphysical pole 
to infinity, but the Kondo regime defined as a — > is 
not properly reproduced. Moreover, FLEX-type approx- 
imations completely smear out the large- frequency struc- 
tures. There are no satellite Hubbard bands j2£ 

A better systematic alternative how to move the RPA 
critical interaction U c to infinity is to introduce a two- 
particle self-consistency. The most straightforward way 
to do so is to use the parquet approach. Its basic 
idea is to replace the bare interaction in Bethe-Salpeter 
equations (|7J| and © by irreducible vertex functions 
A 6 ' 1 and A ee , respectively. In the exact theory then 
$ [A e,i ] = VP [A ee ] = T and equation (JT2J) transforms to 
r + U = A eh + A ee . When we combine the two Bethe- 
Salpeter equations with the parquet equation modi- 
fied as above we obtain a set of non-linear integral equa- 
tions for the irreducible vertices A eh [U, G] , A ee [U, G] and 
finally also for the full vertex T[U,G]. It is, however, 
impossible to solve the parquet equations in the strong- 
coupling regime with a vanishing Kondo scale a — > 0. 
The problem lies in that all the vertex functions, solu- 
tions of the full parquet equations, generally depend on 
three frequencies and there is presently no method, ei- 
ther analytic or numerical, that would predict the full 
analytic structure of the two-particle vertices. We hence 
have to approximate the full parquet equations. 

We proposed a simplification of the parquet equa- 
tions in the Kondo regime of the single-impurity Ander- 
son modeh 21 ' 22 This simplification utilizes a partial sep- 
aration of large and small frequency fluctuations in the 
Kondo regime when the solution of the Bethe-Salpeter 
equation in the electron- hole channel is almost singular. 
In the non-renormalized perturbation expansion it is the 
RPA pole. The principal idea of our simplification is that 
we separate large fluctuations diverging at the critical 
point from those remaining finite. In this way we de- 
scribe correctly universal quantities connected with the 
critical point in the electron-hole Bethe-Salpeter equa- 
tion. We neglect all finite fluctuations and keep only the 
relevant ones diverging at the critical point. We then 
replace all finite frequency-dependent functions by con- 
stants and keep only the frequency dependence in the 
variables controlling large fluctuations. This simplifica- 
tion is kind of a mean-field (simplest) approximation for 
the dynamics in the strong-coupling regime. 

We know from the analysis of multiple electron- 
electron and electron- hole scatterings that only the latter 
can cause a singularity in the vertex function. It means 
that the irreducible vertex in the electron-hole channel 
A eh remains finite. We hence can replace it by an ef- 
fective interaction that we denote U. The rcnormalizcd 



Bethe-Salpeter equation in the electron-hole channel is 
RPA-like with the following replacement 



U 



(14) 



The effective interaction U will be determined from the 
Bethe-Salpeter equation in the electron-electron channel 
by using the parquet self-consistency. Within the scope 
of our approximation the Bethe-Salpeter equation in the 
electron-electron scattering channel holds only in an aver- 
aged form smearing out finite frequency differences. We 
proposed in Ref. [22| a decoupling of the Bethe-Salpeter 
equation in the electron-electron scattering channel so 
that an integral equation is converted to an algebraic, 
solvable one. In this approximate simplification an aver- 
aged integral kernel is assumed to diagonalize the inte- 
gral equation. Such a decoupling is not uniquely defined 
and should be chosen so that the effective interaction is 
a real self-adjoint (positive definite) matrix. An equa- 
tion for the effective interaction derived from a decou- 
pling with the right-averaged integral kernel used in the 
single-orbital case reads 



Uxp = Uip-L* 



1 + L* 



(15) 



where we denoted 



{iv 

—n—n' )G, cri (iui n ')G VC r 2 (iuj- n ') . (16) 



This one-sided averaging in the multi-orbital case can 
lead to a matrix of the effective interaction that is not 
self- adjoint. To avoid this spurious behavior we intro- 
duce also a left-averaged integral kernel in the electron- 
electron Bethe-Salpeter equation. We then define an- 
other equation for the effective interaction 



ipU = i/)U-£* 



l + £* 



C 



(17) 



with a left-averaged integral kernel 



= \ J2 Gact i ( iUJn ' ) G " CT = ) A te (W-n-n' ) 



(18) 



We demand that the matrix of the effective interaction 
U be self-adjoint. It is always the case if we use a sum of 
equations (| 15[) and (|17p to determine the effective inter- 
action. That is, we choose an equation symmetric with 
respect to right and left averaging of the integral kernel of 
the Bethe-Salpeter equation in the nonsingular electron- 
electron channel. 

The dominant contribution to the integral kernel in 
the electron-electron scattering channel comes from the 
singular part of the vertex from the electron-hole scatter- 
ing channel, being A = <!>— U. Since the averaged integral 
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kernels L and C depend on a fermionic Matsubara fre- 
quency, we assume that Eqs. (fT5j). (fT7)) are valid only for 
the lowest frequency, being zero at zero temperature. At 
finite temperatures we use the lowest or a few low fre- 
quencies, symmetric around zero, and add all equations 
for these frequencies. The effective interaction is then 
determined form an average over Eq. (|15|) and Eq. (| 1 7[) 
for a few lowest-lying fermionic frequencies. The low- 
temperature asymptotics is in this way well reproduced. 
The higher the temperature the more Matsubara frequen- 
cies we must take into account to maintain accuracy of 



the zero-temperature solution. We know, however, that 
thermal fluctuation suppress the Kondo resonance and 
smear out the RPA singularity in the electron-hole cor- 
relation function^ 

Equations 0, ([To ]l -([T8 | are our simplified parquet 
equations determining the two-particle vertices $ and U . 
The full two-particle vertex T is determined from equa- 
tion (|12p . Within our approximation it reduces to the 
renormalized RPA vertex <!>. We use it to determine 
the self-energy. The respective Schwinger-Dyson equa- 
tion reads 



m a'v m v 

~ ]p X) X]X U ^ ,X <y' a G L , a {iu n+m )G K , al (iw n ,)Gx', a '{i^n^ m )^^f (i"m) , (19) 

n',m a' i' .k' 



where conservation of the orbital index A' = i + n' — v' 
holds. 

The approximation is almost complete. The last thing 
we have to do is to specify the one-electron propagators 
in the parquet equations for the two-particle vertices. We 
know that the full one-particle self-consistency where the 
one-electron propagators are renormalized by the self- 
energy from Eq. (fl9| does not lead to the Kondo scale 
as U — > oo and also smears out the satellite Hubbard 
bands 32 We hence use the Hartree one-electron propa- 
gators that in the dynamical mean-field theory is 

G ua (iw n ) 

= [ djE ^r rr^ (20) 

J iu> n + [i, + ah - e - ^ J2\,<t> Uu,<f n Xa > 
and for the impurity Anderson model 

G ua (iu> n ) 

1 

iui n + fi L + ah- A(iuj n ) - \ Y^x a' Uu X f' n X a' 

(21) 



We use a general two-orbital interacting Hamiltonian 
comprising also Hund's exchange coupling^ 

Hi = U 2J n- t t n 4 + X [Uir)^ + U 2 S aa >] n lrT n 2 a> 

t go' 
(^l c 2t C H C 2J. C lt ^ 2C 2t C 2i C lJ, C lt H-C.^j (23) 

where T] aa i + 8 aa i = 1. A graphical representation of sin- 



We denoted fj, L = /x — E L , where E L is the center of the 
i's orbital and fi is the chemical potential fixing the total 
charge density. The Hartree propagators are statically 
self-consistent, that is, the particle densities n\ a are de- 
termined from the fully renormalized propagator 

n\ a = 4 X Gxa ( iuJn ~ A£ Ac r(iw„)) . (22) 



We denoted AE the dynamical part of the self-energy, 
a correction to the Hartree term. The chemical poten- 
tial is then adjusted so that the total particle density 
is fixed. The approximate scheme based on the simpli- 
fied parquet equations is now complete. We demonstrate 
explicitly applicability of this approximation on a two- 
orbital model in the next sections. 



IV. TWO-ORBITAL MODEL 



I 

gle terms contributing to this two-orbital Hamiltonian is 
presented in Fig. |6] Each input parameter of this Hamil- 
tonian will be renormalized independently. 

Only one of the bare interacting terms denoted U 2 is 
spin triplet. This term does not mix with the other scat- 
tering processes in the Hamiltonian from Eq. (J23J) . The 
simplified parquet equations for this term are identical 
with those from the single-orbital case. The effective in- 
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FIG. 6. Different elementary scattering processes in the two-orbital model. Each of these processes is renormalized indepen- 
dently. 



teraction C/ 9 is 



U 2 =U 2 



1 



11, a 
J 22,a 



(24) 



with the averaged integral kernel from the Bethe-Salpeter 
equation in the electron-electron channel 



L 



11,(7 



= ^^G la (iuj n ,)G2a{iu~n>)Al 1 2 ]l{iV-n-n>) ■ (25) 
n' 

The dynamical part of the (reducible) vertex from the 
electron-hole scattering channel is 



= —tt4t 



(26) 



The other terms in the interaction Hamiltonian from 
Ecj. (|23|) mix in the parquet equations and must be 
treated simultaneously. Since only specific inter-orbital 
processes are allowed we can simplify the representation 
of the interaction term from the preceding section and 
reduce super-matrices (tensors) with four indices to reg- 
ular matrices with two indices. A relation between ma- 
trix indices and super-indices can be chosen as follows 
1 = (1,1), 2 = (2,2), 3 = (1,2), 4 = (2,1). If we go 
over from super-matrices to matrices we have to distin- 
guish two representations corresponding to two matrix 
multiplications used in constructing the Bethe-Salpeter 
equations. The interaction matrix suitable for summa- 
tion of multiple electron-electron scatterings is 



U ee = 



U J 2 

J 2 U 

Ux Ji 

Ji Ux 



(27) 



The interaction matrix for the electron-hole multiple 
scatterings interchanges the two Hund couplings Ji and 
J2. Hence we have 

fUJx 
Ji U 


V 

Each matrix representation is chosen so that the scatter- 
ing processes in the respective channel are simple matrix 
multiplications. 



U'" : 




(28) 



The matrix integral kernel for the Bethe-Salpeter 
equation in the electron-hole channel is 



K=U X 



( lhxx\ J_ixl \ 

Jix\ tW _o _o 

Uxxl J_2X 2 x 

\ J 2X \ Uxxl J 



(29) 



while the right-averaged kernel in the Bethe-Salpeter 
equation in the electron-electron channel reads 



U 



/ Til rl2 

/ ^11 ^12 

r21 r22 

-^21 n 22 





V 






^22 
7-21 
^12 



\ 



T 12 
-^21 

Lfx) 



(30) 



with 



= \ X] A tv a a2 [W-n-n') G uri (iUn>)G V(72 {iU- n ') (31) 



and A = $ — U, Analogously we construct the left- 
averaged kernel C, if necessary. The matrix form of the 
simplified parquet equations then is in the electron-hole 
channel 



$ = 



1 + K» 



-1 j^_eh 

u 



(32) 



and in the electron-electron channel with right and left- 
averaged integral kernels 



Uip + ipU = Uip- L ee * 



1 + L ee * 



1 + £ ee * 



(33) 



All matrices appearing in the parquet equations (|32[) 
and (|33| are block diagonal. We denote 



$ low 



JJ _ ( C^upp 







(34) 



and solve the parquet equations for each block separately. 
The solution in the electron-hole scattering channel has 
the following representation 
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$upp = 



1 



u i + (u 01 Uo2- 4)xi J i 



i + u aiX \ + u 02X 2 2 + (u 01 u 2 - J') x\xl V Jl u ° 2 + ( c/oiC/ ° 2 _ J ') x} 



(35) 



and 



i ( u 1 + (u\- 4) x \ J 2 

$low = V _ ' _ /_ 2 _ 2 \ | • (36) 

i + c/i (xl + x?) + (pi Jl) xlxi V j2 c/l + I 171 ~ 1/2 ) * 2 

Notice that the Hubbard interaction U splits into two renormalized values t/oi and U02 if the bubbles in different 
orbitals are different, \l xi- Note that we suppressed the frequency index at the electron- hole bubble leading to a 

frequency dependence of vertex $. 

We now choose the singular part of the renormalized electron-hole vertex A = $ — U and use it to obtain a solution 
of the Bethe-Salpeter equation in the electron-electron channel with an averaged (right, left) kernel that determines 
the matrix of the effective (renormalized) interaction. Its upper matrix block reads 



r- , EWl J 2^2 \_( Ui>\ J 2 ^l \ 1 

UuppV— -j ,1 77 ,2 - t -1 TT i2 



Jii>\ U m i>i J \ J 2 ^l U^l J (1 + Lff )(1 + - Lf§L§f 

ree ree /i r ee \ j ee ( j ee \ j ee \ j ee t ee t ee j ee\ 

^12^21 I 1 ^li; ^12 1^11 ^22 T" ^11^22 — ^12^2lJ 

t ee t ee j ee t ee\ /i i j ee \ I j ee \2 , ree ree /i j ee \ 

^11^22 ~" ^12^217 I 1 "T MlJ 1^22; ~T~ ^21-^12 I 1 ~~ -^227 




(37) 



The lower block analogously is 

fTli, = ( unit Jirt \_{ 

UlowiP- -J , 1 77 ,2 - j ,1 T r, 2 



j^i Uirtj v j i^2 vM) (1 + + lid - mLf 3 

(/1 1 r ee \ / r ee \2 i ree ree /i ree\ j ee ( j ee \ j ee \ j ee j ee j ee j ee\ 

I 1 ' -^44J 1^337 ^34^43 I 1 ~~ ^33^ ^34 1^33 ^44 "T ^33^44 — ^34^43^ 

T ee f r ee i ree i t ee t ee j ee j ee\ /i i j ee \ ( j ee \* , ree ree /i ree\ 

^43 1^33 "+" -^44 ^33^44 ~~ ^34^437 I 1 + ^33 J 1^447 + ^43^34 I 1 — ^44J 



(38) 



The electron-electron bubble ip and vertices L ee are taken 
here (zero temperature) only at zero frequency. If ip\ ^ 
ip2 we must symmetrize Eqs. ([37]) and (|38|) by adding the 
Bethe-Salpeter equation with the left averaged kernel C 
so that to keep the matrix of the effective interaction 
self-adjoint. 

We use the electron-hole dynamical vertex $ result- 
ing from the simplified parquet equations and calcu- 
late the self-energy correction to the Hartree term from 
Schwinger-Dyson equation (fn?|) . It is in the super-matrix 
representation 

AE La (iui n ) = - - u t\° a <xV^\'Z° i( i v ™) 

m A, A' l',er' 

x G\a>(iu} n + iv m ) . (39) 

The dynamical part of the self-energy AE is then used 
to determine the particle densities from Eq. ([22)) so that 
the static one-electron self-consistency is fulfilled. 



V. RESULTS 

We analyze a general two-band Hubbard model within 
the dynamical mean-field theory. We set the model pa- 
rameters U\ = U — 2 J, U2 — U — 3 J and J\ = J 2 = J, 
which corresponds to a spherically symmetric situation. 38 
We choose the Hund exchange so that all three electron 
interactions remain positive. We use the one-electron 
propagator from Eq. (|2Tj)) with the semi-elliptic density 
of states, pi(e) = p 2 (e) — 2k~ x \J\ — e 2 the bandwidth 
of which is w = 2. We investigate only the spin sym- 
metric case where the one-electron propagators are spin- 
independent. The critical interaction in the random 
phase approximation, separating the weak and strong 
coupling regimes, then is U c = 3%/8 w 1.18. It does not 
depend on the number of orbitals if the Hund coupling 
J = 0. We perform all calculations at zero temperature. 

We first investigate the simplest situation of the 
charge-symmetric model at half filling, that is n = 2, 
with no orbital splitting A = \E% — E\\/2 = and with 
vanishing Hund coupling, J = 0. This situation mimics 
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the single-orbital model. The emergence of the Kondo 
scale is plotted in Fig. [7J The quasiparticle peak starts 
to develop from the bare density of states for interac- 
tions around the critical one from RPA, U c « 1.18. The 
satellite Hubbard bands are well formed in the strong- 
coupling regime U > U c and lie slightly beyond the 
atomic values ±{7/2. The Kondo scale remains expo- 
nentially small but non-zero even for large interaction 
strengths. The Kondo asymptotics for the semi-elliptic 
density of states is a a exp{— AUpa/3}, where po is the 
density of states at the Fermi energy. This result is in 
discrepancy with general confidence that the Hubbard 
model at half-filling undergoes in the paramagnetic phase 
a Mott-Hubbard metal insulator transition at a finite in- 
teraction strength. The existence of a metal-insulator 
transition in DMFT is concluded from extrapolations of 
low-temperature Monte- Carlo£ or NRG^^ data. Al- 
though there is no doubt about the existence of an insu- 
lating solution for very large interaction strengths match- 
ing the atomic limit, numerical methods cannot presently 
exclude the existence of a metallic solution with a very 
narrow (below numerical resolution) quasiparticle peak. 
The existence of a sharp metal-insulator transition has 
not yet been rigorously proved and there is at present 
neither a consistent solution at the critical point of a 
transition nor there is an analytic formula for the vanish- 
ing width of the quasiparticle peak. If there is a metal- 
insulator transition, the metallic side of the transition 
cannot be a Fermi liquid. The low-frequency asymp- 
totics of the electron-hole correlation function in Fermi 
liquid is due to sum rules of order uj^ 1 . Such a behavior 
is then incompatible with Bethe-Salpeter equations that 
have to be fulfilled self-consistently in the electron-hole 
and electron-electron scattering channels by an exact so- 
lution. Only integrable singularities can emerge in an 
exact solution and Fermi liquid hence does not allow for 
a divergence in the electron-hole correlation function ex- 
pected at the Mott-Hubbard metal-insulator transition. 
In this respect the simplified parquet approximation is 
consistent with the exact two-particle behavior. 

Nonexistence of a sharp metal-insulator transition in 
the simplified parquet equations is not caused by ap- 
proximations made in this approach but rather it is 
a consequence of a Fermi-liquid solution fulfilling self- 
consistently Bethe-Salpeter equations simultaneously in 
the electron-hole and electron-electron scattering chan- 
nels. There is an insulating paramagnetic state in 
parquet-based theories for U > U c (critical interaction 
of the RPA pole) being a superposition of Hartree po- 
larized solutions, but it coexists with the metallic one in 
the strong-coupling regime and has a higher total energy. 
A metal-insulator transition in the Fermi liquid can ei- 
ther be connected with a breaking of a global symmetry 
and/or can be noncritical as we discuss later for orbital 
splitting in a crystal field. 

We compared the spectral function calculated from 
the simplified parquet equations with other diagram- 
matic approximations based on summations of Feynman 
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FIG. 7. (Color online) Emergence of the Kondo scale in the 
two-orbital, charge, spin and and orbital symmetric model 
without Hund exchange (one electron per orbital). Semi- 
elliptic density of states. 



diagrams for multiple two-particle scatterings. We chose 
a value of the interaction strength U = 3 so that to be 
in the strong-coupling regime where renormalizations of 
the perturbation expansion are necessary. We used the 
spectral function for the self-consistent T-matrix, iter- 
ated perturbation theory in second order (IPT) and for 
T-matrix with IPT (topological) self-consistency, that is, 
the local renormalizcd propagator of the Hubbard model 
G is defined as Q' 1 ^) = G _1 (w — E(w)) +E(w). The re- 
sults of the comparison are plotted in Fig. [5J We can see 
that except for the fully self-consistent T-matrix, the cho- 
sen three approximations display the expected three-peak 
structure. They differ, however, in the way the satel- 
lite bands are attached to the central quasiparticle peak. 
The simplified parquet approximation keeps a small finite 
density of states between the central and satellite bands. 
It means that this approximation allows for an energy 
transfer between the central and the satellite bands. Ap- 
proximations based on the IPT self-consistency develop 
a quasi-gap between the low and high energy states and 
no energy exchange occurs. The satellite bands are less 
spread in the parquet approximation than in the IPT- 
based ones. The detailed shape of the central quasi- 
particle peak in these approximations is magnified in 
Fig. IH1 The most narrow is the central peak for IPT 
suggesting vicinity of a metal-insulator transition with 
zero width of the quasiparticle peak. The parquet-based 
approximations, on the other hand, slow down the ten- 
dency toward a metal-insulator transition. They actually 
screen it, since there may be no critical metal-insulator 
transition in these theories as discussed above. 

The simplified parquet equations show best agreement 
with the exact solution of the single-impurity Anderson 
model at half filling, where the density of states at the 
Fermi surface is fixed. The same holds also for other 
fillings at intermediate and strong coupling calculated 
within the dynamical mean-field approximation. Apart 



10 



U=3 



U=2 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




1 1 1 


parquet 

TMA(IPT) 

SCTMA 


i 
i 

/ \ I 


\ IPT 

i / \ 


1 \ ' 

/ A ' 

.....J S\ ' 
!**• (' 




\* V = VN - ; 

,1 : \ v 


\ //jv 'A 


'I " \ / 

• ■ i V.-'\ J 




..-U ' \ ; r'i 4- 



-2 



-1 1 
Energy 



O 
Q 




FIG. 8. (Color online) Spectral function for the symmetric 
situation calculated in various diagrammatic approximations: 
simplified parquets, self-consistent T-matrix (SCTMA), iter- 
ated perturbation theory (IPT), and T-matrix with IPT prop- 
agators. 
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FIG. 9. (Color online) Detailed structure of the quasiparticle 
peak for the approximations from Fig. [8] 



from the charge-symmetric situation the simplified par- 
quet equations reproduce qualitatively well the low and 
high frequency features of the spectral function, although 
the width of the quasiparticle peak and the Kondo tem- 
perature derived from it do not follow as closely the ex- 
act solution for the flat-band model as for the half-filled 
case<^ The spectral function of a strong-coupling solu- 
tion changes when receding from half filling so that the 
central quasiparticle peak broadens and slightly moves 
below the Fermi surface. The lower satellite band moves 
toward the Fermi energy. The lower band and the central 
peak eventually merge, the quasiparticle peak is absorbed 
by the satellite one, when the electron-hole asymmetry 
becomes prominent. This feature seems to be universal 
for all theories and is demonstrated in our approach in 
Fig. [TO] 

We hitherto have shown the results for the two-orbital 



FIG. 10. (Color online) Spectral function in charge asym- 
metric situations for various filings n, averaged number of 
electrons per orbital. 



model without Hund exchange and orbital splitting, that 
is for J = and A = 0. Such a situation is close to 
the single-orbital model as we demonstrated. The most 
important change in the model with nonzero Hund cou- 
pling is a splitting of the interaction strength, that is, 
U\ = U — 2J and U2 = U — 3J. It means that the 
interaction strength controlling the approach to the crit- 
ical point in the electron-hole channel is the Hubbard U 
between different spins on the same orbital. The Hund 
exchange decreases the critical interaction at which the 
electron-hole vertex diverges in the bare perturbation 
theory. It would then narrow the width of the Kondo res- 
onant peak at a fixed renormalized Coulomb repulsion U. 
In the real situation we, however, do not fix the renormal- 
ized but rather the bare Coulomb repulsion. We find that 
though larger Hund exchanges lead to higher values of the 
averaged kernels L ee , as we can see from Eqs. (|37|) and 
they simultaneously cause decrease in the renor- 
malized repulsions U a at a fixed bare Coulomb interac- 
tion. Consequently, the weight of the central quasipar- 
ticle peak in the self-energy decreases. We plotted the 
impact of the Hund coupling on the spectral function in 
Fig. QTJ We observe that the width of the Kondo peak 
shrinks but only on a narrower frequency interval near 
the Fermi energy. The Hund coupling at a fixed bare 
Coulomb repulsion then leads to filling of the states be- 
tween the central peak and the satellite bands. Further, 
the satellite bands are positioned closer to the central one 
and become less pronounced with the increasing Hund 
exchange. 

The most interesting on the two-orbital model are 
changes in the spectral function due to a broken or- 
bital symmetry induced by a crystal field. That is, when 
the one-electron propagators are different on different or- 
bitals. We achieve this here by splitting the centers of the 
two bands. We set E\ = — E 2 = A, but keep the same 
band structure, pi(e) — /92(e)- The effect of the orbital 
splitting depends on whether we have a finite or infinite 
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FIG. 11. (Color online) Effect of the Hund exchange coupling 
on the spectral function. 

bandwidth of macroscopically occupied energies. Orbital 
splitting leads in the former case to a metal-insulator 
transition unlike the latter one. The genesis of a metal- 
insulator transition in the model with the semi-elliptic 
density of states is plotted in Fig. Q21 We chose a mod- 
erate interaction strength, U = 1.5, so that to be able 
to follow the process more closely. When we increase 
the orbital splitting parameter A there is almost no ob- 
servable change at the Fermi energy with a pronounced 
quasiparticle peak. It slightly splits into two but stays at 
the Fermi energy. More apparent is the change caused 
by orbital splitting in the bulk of the energy band where 
the states are rearranged so that the orbital polarization 
is enhanced. Such a response on the band splitting ef- 
fectively pushes the system out from the critical region 
of the RPA pole. With the increasing band splitting the 
dynamical renormalization of the band structure, apart 
from the vicinity of the Fermi energy, becomes less impor- 
tant and the solution approaches the Hartree one. The 
process continues up to the point when one of the or- 
bital saturates and the other is emptied. The insulating 
solution is then that of the Hartree approximation with 
no dynamical correlations. Changes in the quasiparticle 
peak, split and a broadening, are observed only close to 
the transition. They vanish practically at the transition 
point. There are, however, no singular two-particle func- 
tions at the transition point, the metal-insulator tran- 
sition is hence noncritical. It is worth noting that the 
Kondo quasiparticle peaks exist in a pronounced way 
only in the strong-coupling regime, cf. Fig. 1131 

We further studied the effect of the Hund exchange on 
the orbital metal-insulator transition. We set the param- 
eters of the model so that we could compare our results 
with the existing Monte-Carlo simulations of Ref. H3- 
We chose three values of the bare orbital splitting 2A = 
0.2, 0.4, 1.0. For each value of the band splitting we fixed 
four ratios of the Hund coupling J/U = 0,0.05,0.1,0.2 
and plotted the filling (per spin) of the upper orbital 
m as a function of the bare interaction at half filling 




-3-2-10123-3-2-10123 



FIG. 12. (Color online) Orbital splitting and a metal- 
insulator transition in the half-filled case for the semi-elliptic 
density of states and U = 1.5. 
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FIG. 13. (Color online) Orbital splitting and a metal- 
insulator transition in weak coupling, U = 0.5. 



ni + 7i2 = 1 in Fig. [TU We observe that small Hund 
couplings do not prevent the system from falling into the 
orbitally polarized state as demonstrated in Fig. [12j For 
large A the occupation of the upper orbital decreases 
practically monotonically with the increasing interac- 
tion. For smaller splittings the filling n\ increases and 
may get to the Kondo regime at intermediate interaction 
strengths where both orbitals are almost equally occu- 
pied (ni w ri2 ~ 0.5). Unlike the Monte-Carlo data, we 
do not reach, however, the Mott insulating phase from 
the Fermi-liquid solution. If J/U < 0.2 the system does 
not stay in the Kondo regime for ever. Strong electron re- 
pulsion increases orbital polarization due to the Hartree 
shift and the system eventually goes over into a polar- 
ized insulator n% = 0. For ratios J/U > 0.2 the Hartree 
term starts to act against the crystal field and diminishes 
orbital polarization. The filling m then asymptotically 
approaches half filling. Apart from nonexistence of the 
Mott transition these effects of the Hund coupling on 
the occupation of orbitals in a crystal field are described 
by the simplified parquet equations in a good agreement 
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FIG. 14. (Color online) Filling per spin of the upper orbital 
as a function of the bare interaction strength U. Different 
curves for a given split 2A correspond to different J values 
from bottom to top J/U = 0, 0.05, 0.1, 0.2. A metal-insulator 
transition (ni = 0) takes place for values of J/U < 0.2. The 
Kondo regime, vicinity of the Mott-Hubbard transition (ni ~ 
1/2) is reached only for larger values of J/U. 
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FIG. 15. (Color online) Orbital splitting in the single- 
impurity Anderson model for U — 5. No metal-insulator 
transition takes place. 



with the Monte-Carlo simulations. It is worth mentioning 
that the Kondo resonance for J/U < 0.2 survives to a cer- 
tain value of the interaction strength beyond which it is 
rather abruptly suppressed, similarly to the scenario from 
Fig. [121 Vanishing of the Kondo resonance is followed by 
a sharp saturation of orbital polarization. The numerical 
solution then becomes unstable, since the metallic phase 
in the extreme Kondo regime and the orbitally polarized 
insulating solution coexist with almost the same energies. 

The scenario of orbital splitting changes if the bare 
energy bands are infinite. When an orbital symmetry 
breaking is switched on, the system with infinite band- 
widths behaves for small splittings analogously to the 
model with finite energy bands. The changes take place 
in the bulk, away from the Fermi energy so that the or- 



bital splitting is enhanced, cf. Fig. [T5] Similarly there 
are almost no observable changes on the quasiparticle 
peak. The system is slowly pushed away from the criti- 
cal region of the RPA transition by increasing the orbital 
splitting and approaches orbital saturation. But due to 
the infinite band width, there is no transition to the 
Hartree insulator and the density of states at the Fermi 
energy remains nonzero for all interaction strengths. The 
Kondo peaks do not survive to infinite interaction. A 
sudden macroscopic cleft breaks in between the central 
quasiparticle peaks of different orbitals at a finite inter- 
action strength. The split Kondo peaks are then soon 
absorbed by the bulk energy bands. The solution ap- 
proaches a weak-coupling (Hartree-Fock) state with al- 
most separated energy bands centered around energies 
±(A + Um /2), where m is the orbital moment. 



VI. DISCUSSION AND CONCLUSIONS 

Dynamical effects of strong electron correlations in 
impurity models manifest themselves in a three-peak 
structure of the spectral function with a narrow central 
quasiparticle peak near the Fermi energy and two satel- 
lite bands centered around eigenvalues of the local part 
of the interacting Hamiltonian. There is no exact solu- 
tion proving such a picture, but a number of approximate 
numerical and analytic-numerical solutions confirm this 
behavior of the spectral function. There are a number of 
such approximate schemes for the simplest single-orbital 
impurity and DMFT models. Extensions of approxima- 
tions showing the expected three-peak spectral function 
in strong coupling to multi-orbital models face a consider- 
able increase of complexity in the strong-coupling regime. 
Monte-Carlo simulations as well as the numerical renor- 
malization group applied on multi-orbital models need 
to resort to additional approximations to reach sufficient 
accuracy of dynamical and spectral properties. 

Most of the existing approximations on dynamical 
properties of strongly-correlated electron systems make 
directly accessible only one-electron properties. Only 
few of them, even in single-orbital models, are able to 
approach two-particle quantities and in particular, to 
control two-particle divergencies that develop in Bethe- 
Salpeter equations. We extended in this paper an ap- 
proximation based on parquet equations controlling the 
RPA singularity in the electron-hole scattering channel 
to a general multi-orbital model of correlated electrons. 
Approximations used to simplify complexity of the full 
parquet approach resulting from a separation of large and 
small dynamical fluctuations in the critical region of the 
RPA pole allow us to reproduce the expected three peak 
structure of the spectral function of SIAM in the Kondo 
regime. An extension of this approximation from single 
to multi-orbital models proved to be manageable with 
only a moderate increase in numerical requirements. An 
Z-orbital model is in the most general case described by 
4l 2 x 4i 2 matrices and the numerical expenses increase 



13 



polynomially with the number of orbitals. In real situ- 
ations symmetries and conservation laws further reduce 
complexity of the simplified parquet equations. 

The impurity solver for the dynamical mean-field the- 
ory we built up within the parquet approximation suits 
best for describing the dynamical and spectral properties 
of the metallic (Fermi liquid) phase from weak up to mod- 
erately strong electron correlations. No self-consistency 
is needed in the weak-coupling regime below the criti- 
cal interaction strength U c at which the RPA vertex di- 
verges. The parquet self-consistency is indispensable in 
intermediate (U s» U c ) and strong (U > U c ) coupling. 
The numerical procedure for the solution of the simpli- 
fied parquet equations gets unstable for large interaction 
strengths. We were able to reach a stable metallic solu- 
tion up to U ~ 5U C . Unlike other approaches such NRG, 
QMC or IPT this approximation does not allow for a 
critical metal-insulator transition without a symmetry 
breaking in Fermi liquid. This is a general feature of 
all impurity solutions where Bethe-Salpeter equations in 
the electron-hole and electron-electron channels should 
be obeyed simultaneously in a self-consistent manner. 
Parquet equations are hence unable to describe a Mott 
insulator without a symmetry breaking. The insulating 
solution obtained within the local parquet approach is a 
superposition of fluctuations-free Hartrcc saturated solu- 
tions and the transition to such a state is noncritical. 

The simplified parquet equations presented in this pa- 
per are not quantitatively as accurate as numerically ex- 
act solutions, but they offer a numerically manageable 
scheme with an analytic control of two-particle singular- 
ities in Bethe-Salpeter equations. The self-energy and 
the spectral function calculated within this approxima- 
tion obey Fermi-liquid rules in the metallic phase and 
produce the typical three-peak structure in the density 
of states in the strong-coupling regime. The width of 
the central quasiparticle peak is exponentially small for 
large interaction strengths. This approximation may not 
reproduce the exact asymptotics of the Kondo scale and 
temperature in the single-impurity model as well as e. g. 
the local moment approach, but unlike the latter it does 
not show any unphysical behavior and does not use auxil- 
iary fitting parameters. The approximation based on the 
simplified parquet equations is formulated in such a way 
that it can directly be applied not only to impurity mod- 
els but to lattice models within the dynamical mean-field 
approach as an impurity solver. 

The impurity solver based on the simplified parquet 
equations is a semi-analytic theory for local spectral 
properties formulated in real frequencies with a full con- 
trol of analytic properties of the approximation. From 
this reason the approximation best suits to the calcula- 
tion of spectral properties and the dynamical self-energy 
directly accessible from a perturbation theory. Its unique 
feature is that it offers analytic control of two-particle 
singularities in the Bethe-Salpeter equations. The sim- 
plified parquet approximation is well suited to problems 
involving two-particle dynamical, and in particular criti- 



cal, functions. The approximation is less suited to global 
thermodynamic quantities being averages of dynamical 
functions over frequencies. An extension to finite tem- 
peratures does not represent a problem and the effective 
interactions can be calculated also at nonzero tempera- 
tures as discussed in this paper. The quantitative relia- 
bility of the simplified parquet equations decreases with 
the reduced influence of the critical behavior of the vertex 
function in the vicinity of the RPA pole (Kondo regime) . 
Only in the critical region of a pole in the two-particle 
vertex we can rely on a separation of large and small 
dynamical fluctuations, that is, a decoupling of singular 
and non-singular functions. 

To our knowledge, the simplified parquet equations is 
the only scheme that uses a two-particle self-consistency 
to control the Kondo asymptotics and criticality of the 
RPA pole. There are attempts to use parquet equa- 
tions for renormalizations of the perturbation expansion, 
but they are based on the full one- and two-particle 
self-consistency 4 1 1 42 It is, however, important to realize 
that the full dynamical one-particle self-consistency im- 
pedes control of two-particle divergencies and the Kondo 
scale in the spectral function is no longer correctly repro- 
duced. Moreover, the full dynamical one-electron self- 
consistency ultimately leads to smearing of the satellite 
bandsi^i It is hence of importance that the parquet equa- 
tions are solved with the Hartree propagators so that to 
reproduce the three-peak spectral function with the cor- 
rectly scaled central quasiparticle peak. 

The extension of the simplified parquet equations was 
used in this paper as an impurity solver for the two- 
band Hubbard model with Hund coupling and orbital 
splitting. We showed emergence of the Kondo quasi- 
particle peak with the increasing interaction strength at 
half filling. The quasiparticle peak survives only close 
to the electron-hole symmetry where multiple electron- 
hole scatterings are dominant. When deviations from 
this symmetry are significant the quasiparticle peak is 
absorbed by one of the satellite bands. We also demon- 
strated how a metal-insulator transition induced by a 
crystal field takes place in a model with finite energy 
bands. The parquet equations do not allow for a criti- 
cal transition in the local Fermi liquid. We showed that 
the orbitally induced metal-insulator transition is non- 
critical, but with a Kondo behavior almost up to the 
transition point, if the interaction is sufficiently strong. 

Last but not least, the approximate impurity solver 
we built up on parquet equations uses the formalism of 
Green functions and can conveniently be combined with 
ab-initio computational schemes to assess correlation ef- 
fects in materials with a nontrivial atomic structure. 
It should be suitable in particular for metallic systems 
where we expect Kondo-likc behavior or strong valence 
fluctuations. 
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